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ABSTRACT 

BL Lacertae is the prototype of the BL Lac class of active galactic nuclei, exhibiting 
intensive activity on parsec (pc) scales, such as intense core variability and multiple 
ejections of jet components. In particular, in previous works the existence of preces- 
sion motions in the pc-scale jet of BL Lacertae has been suggested. In this work we 
revisit this issue, investigating temporal changes of the observed right ascension and 
declination offsets of the jet knots in terms of our relativistic jet-precession model. The 
seven free parameters of our precession model were optimized via a heuristic cross- 
entropy method, comparing the projected precession helix with the positions of the jet 
components on the plane of the sky and imposing constraints on their maximum and 
minimum superluminal velocities. Our optimized best model is compatible with a jet 
having a bulk velocity of 0.9824c, which is precessing with a period of about 12.1 yr 
in the observers reference frame and changing its orientation in relation to the line of 
sight between 4° and 5°, approximately. Assuming that the jet precession has its origin 
in a supermassive binary black hole system, we show that the 2.3-yr periodic variation 
in the structural position angle of the very-long-baseline interferometry (VLBI) core 
of BL Lacertae reported by Stirling et al. is compatible with a nutation phenomenon 
if the secondary black hole has a mass higher than about six times that of the primary 
black hole. 

Key words: galaxies: active - BL Lacertae objects: individual: (BL Lacertae) - 
galaxies: jets - radio continuum: galaxies. 
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1 INTRODUCTION 



BL Lacertae or 2200+420 (z=0.0686; IVermeulen et all 
1 19951 ). the prototype of the BL Lac class of active galactic 
nuclei, is hosted by an elliptical galaxy, composed mainly by 
a stellar population of about 0.7 Gyr (Hyvoncn et al. 20q3)- 
As usual in BL Lac type objects, its nuclear region exhibits 
strong continuum variability on different time-scales, from 
days to years, over the whole electromagnetic spectrum. 

High-resolution interferometric images at radio wave- 
lengths show the presence of a comp act core and a di ffuse 
halo- like source at arc-second scales l|Antonuccilll986l ). At 
an intermediate resolution (~ 10 mas), BL Lacertae shows 
a core-jet structure, with compo nents following differe nt tra- 
jectories on the plane of the sky (iPol atidis et al.iri995l ). Sim- 
ilar behavior is seen at mas resolution, wh ich has been at- 
tributed to either helical instabilitie s (e.g.. pRitevama et al.l 
ll998l:lDenn. Mutel fc Marscheij lioOO) or jet inlet precession 
l|Stirling et al.ll2003l : iTatev amal 12009.1 . Concerning jet pre- 
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cession, IStirling et all l|2003l ) proposed a period of 2.29 years 
based on the periodic variations in the polarization position 
angle at 1 mm, and in the direction of the innermost ra- 
dio core component at 43 GHz, whi ch has been matter o f 
criticism in som e recent works (e.g., iMutel fc DennI [20051 1 . 
iTatevarnal l|2009l ) proposed an alternative precession model 
in which the BL Lacertae parsec-scale jet changes its ori- 
entation in a period of about 26 yr. Such claim was based 
on analyses of maps of BL Lacertae at 8 and 15 GHz in a 
super-resolution mode. 

In this work, we reanalyse the jet-precession proposition 
in terms of our ballistic jet precession model, described in 
Section 2. In Section 3, we introduce the cross-entropy (CE) 
global optimization technique in the context of our preces- 
sion model. In Section 4, we show the general results from 
the application of our CE jet-precession model to BL Lac- 
ertae, as well as the observational constraints used in this 
work. We explore in Section 5 the possible consequences of 
the underlying jet intensity on the radio and optical histor- 
ical light curves. In Section 6 we study the viability of a 
supermassive binary black hole system in the nuclear region 
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of BL Lacertae in producing the inferred jet precession rate, 
as well as a nutation motion with period of 2.3 yr. Finally, 
general conclusions are presented in Section 7. 

We will assume throughout this work a ACDM cos- 
mology with Ho = 71 km s~^ Mpc"\ f^M = 0.27, and 
f^A = 0.73, which implies 1 mas = 1.296 pc and 1 mas yr^^ 
= 4.516c for BL Lacertae. 



2 BALLISTIC JET PRECESSION MODEL 

Let us consider a relativistic jet receding from the core 
with a constant bulk velocity /3 (in units of the speed of 
light c) that precesses around a fix ed axis (see Fig. 1 in 
ICaproni. Monteiro fc Abrahamll2009l for a schematic repre- 
sentation of this). Precession makes the jet inlet direction 
vary with time with a precession period Pprcc.s measured 
in the source reference frame, producing a cone with semi- 
aperture angle ipo- The precession phase uJsAts = 2-K{ts — 
io,s)/Pproc,s = 27r(rs — ro,s) i s chosen arbitrarily to be zero on 
the y s-gs-plane at Ts = ro,s (|Caproni. Monteiro fc AbrahamI 
l2009f ). From these definitions, the instantaneous jet inlet di- 
rection is given in terms of a unit vector with Cartesian 
components: 

ex,s(Ts) = sin sin[t27r(rs - ro,s)], 
ey,s(Ts) — siniy3o cos[t27r(rs — to,s)], 
ez,B{rB) = cos ((30, 

where t gives the sense of precession, t = 1 for clockwise 
sense and l — ~1 for counterclockwise precession. 

Introducing the parameters <j)o, the angle between the 
precession cone axis and the line of sight, and rjo, the posi- 
tion angle of the axis on the plane of the sky (p ositive from 
north to east), we have (e.g., Caproni fc Abrahamii2004a .b: 
ICaproni. Monteiro fc AbrahamI I2OO9I '): 



ex,obs(Ts) = A(ts) cos 770 — ey,s(Ts) sin 770, 
ey,obs(Ts) = ^(ts) sin 770 + ey,s(rs) cosryo, 
ez,obs(Ts) = -ex,s(Ts) sin^o -I- ez,s(Ts) COS0O, 
where: 

A{ts) = ex,s(-rs) cos(/>o -I- ez,s(-rs) sin<;/>o- 



(1) 
(2) 
(3) 

(4) 



The instantaneous angle between the jet and the line of 
sight cf) is calculated from: 



(ts) = arccos[ez,obs(Ts 



(5) 



while the position angle of the jet on the plane of the sky rj 
is obtained from: 



77(ts) = arctan 



,obs V '^s 



6x,obs ('7~s) 

The observed jet velocity /3obs(Ts) is: 
/3obs(-rs) = 7/3(5(rs)sin(^(rs), 
where the jet Lorentz factor 7 is 

and the jet Doppler factor 5 is: 
^(•rs) =7"Ml-/9cos<^(r,)]-\ 



(6) 

(7) 
(8) 
(9) 



In order to compare predictions from the preces- 
sion model with observational data, it is necessary to 
transform the elapsed time measured in the source's 
reference frame dts to the tim e interval in the ob- 
server's framework dtoha- Following Gower et al.l (|l982l ') and 
ICaproni. Monteiro fc AbrahamI (|2009l ). we can write: 



Pprcc.obs S^^[T)dT 

where Atohs = (tobs - io.obs^ 
The relation between 
source's framework and that measured by the observer 



(10) 

and Atb = Ts — To,B- 
the precession period in the 



,,obs IS given as: 



prcCjObs 



(1 + ^) /g'5-i(r)dr' 



(11) 



where z is the redshift of the source. 

A fiuid element of the jet, ejected at time tobs.oj from a 
jet inlet region that is precessing according to our model, will 
be observed at time fobs with right ascension and declination 
offsets from the core Aaniod and A5niod respectively, given 
by: 

AQmod(iobs) = {tabs — iobs.oj ) M(*obs,Gj ) slu [j7(tobs,cj )] , (12) 
A(5inod(iobs) — (tobs — tobs.cj) ^(iobs.cj) cos[r;(tobs,Gj)l, (13) 

where fi is the proper motion of the fiuid element predicted 
by our jet-precession model assuming outward motion at a 
constant speed. 

Considering all the possible values of fobs.cj < fobs we 
obtain a snapshot of the model jet at time fobs, projected 
on the plane of the sky. The real jet, as observed with VLBI 
techniques, is not continuous but formed by discrete compo- 
nents, ejected at unknown epochs with unknown velocities. 
To avoid incorrect identifications of the same component in 
maps obtained at different epochs (which could led to incor- 
rect determination of fobs.cj and /i(fobs,cj) for this compo- 
nent), we will use only the projected position of the compo- 
nents in the plane of the sky for each epoch in which a map 
is available, and compare them with the positions predicted 
by the model. 

In summary, our jet precession model has seven free 
parameters (Pprec,obs, t, 7, Vo, 4>o, fo and ro,s), which are 
determined, via the CE technique as will be discussed in the 
next Section. 



3 THE CROSS-ENTROPY PRECESSION 
METHOD 

3.1 General overview 

The CE method was originally employed in the optimiza- 
tion of complex co mputer simulation models involving rare 
events simulations |^b^stei3[l993), having been modi- 
fied by iRubinsteinI (| 19991 ) to deal with continuous multi- 
extremal and discrete combinatorial optimization prob- 
lems. Its theoretical asymptotic c onvergence in su ch situ- 
ations has been demonstrated by MargolinI (|2004l ). while 
iKroese. Porotskv fc RubinsteiiJ (|20od ) studied the efficiency 
of the CE method in solving continuous multi-extremal op- 
timization problems. Some examples of robustness of the 
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CE m ethod in several situations are listed in Ide Boer et all 
l|2005h . 

CE optimization involves basically random generation 
of the initial parameter sample (obeying some predefined 
criteria) and selection of the best samples based on some 
mathematical criterion. Subsequent random generation of 
updated parameter samples from the previous best can- 
didates are performed iteration by iteration until a pre- 
specified stopping criterion is fulfilled. 

Application of the CE method in astrophysical con- 



texts c an be found i n ICaproni. Monteiro fc Ab: 
tOO^). iMonteiro et al.j (|2010l ). [Caoroni et al.l (|2011 



iMonteiro fc Diaj (|201ll ) 



ahami 



and 



3.2 Our CE algorithm for jet precession modelling 

FoUowing lCaproni. Monteiro fc AbrahamI (12009'). let us sup- 
pose that we wish to study a set of A^d observational data in 
terms of an analytical model characterized by Np parameters 
Pi,P2, ■••iPiVp- In the case of jet precession modelling, the 
observational data correspond to right ascension and dec- 
lination offsets, Aa and AS respectively, for the Ad knots. 
As mentioned before, our precession model is defined by the 
free parameters P (or 7), tjo, 0o, fo, to,s, -Pprec.obs (or Pprcc.s) 
and i, i.e., Np = 7. In this work, we adopted the same pro - 
cedure suggested bv lCaproni. Monteiro fc AbrahamI (|2009l ). 
in which for a given (fixed) value of Pprcc.obs and t, the re- 
maining five parameters are CE optimized. Therefore, it is 
necessary to test different values for the precession period, 
as well as distinct senses of precession in order to obtain the 
best set of model parameters. 

The main goal of the CE continuous multi-extremal 
optimization method is to find the set of parame- 
ters X* — {pi,P2, ■■■,PNp) for which the model pro- 
vides the best description of the d ata (iRubinsteinI 1 19991 : 
iKroese. Porotskv fc Rubinsteinll2006l ). It is performed gen- 
erating randomly N independent sets of model parameters 
X = (xi,X2, ...,xjv), where Xi = {pii,p2i, ...,pNpi), and min- 
imizing a merit function S'(x) used to transmit the quality 
of the fit during the run process. If the convergence to the 
exact solution is achieved then ^(x*) — >■ 0. 

In order to find the optimal solution from CE optimiza- 
tion, we start by defining the parameter range in which the 
algorithm will search for the best candidates: p™'" ^ Pj{k) ^ 
pf^^ , where k represents the iteration number. Introducing 
Pi(0) = (pf " -h pr'')/2 and crj (O) = (pj"'''' - pf ")/2, we 
can compute X(0) from: 



A.,(0) =p,(0) + f7,(0)G.„ 



(14) 



where dj is an A x Ap matrix with random numbers gen- 
erated from a zero-mean normal distribution with standard 
deviation of unity. 

The next step is to calculate Si{G) for each set of Xi(0), 
ordering them according to increasing values of Si. Then 
the first Aeiite set of parameters is selected, i.e. the Aouto- 
samples with lowest S values, which will be labelled as the 
elite sample array X''''*°(0). 

We then determine the mean and standard deviation of 
the elite sample, Pj''*°(0) and crj''''^(0) respectively, as: 



—elite ( r\\ 

P3 (0) 



CTj (0) 



Ae, 



1=1 



(15) 



V [A<=f^(0)-pf-(0)]'. (16) 



^ (Aeiite - 1) ^ 

The array X at the next iteration is determined as: 
^..(l)=pf"(0) + af'^(0)G.,, (17) 

This process is repeated from equation (14), with dj 
regenerated at each iteration. The optimization stops when 
the maximum number of iterations fcmax is reached. 

In order to prevent convergence to a sub-optimal 
solution due to the intrinsic rapid converg e nce o f the 
CE method, |k roese. Porotskv fc RubinsteinI (|2006l ) sug- 
gested the implementation of a fixed smoothing scheme for 

— elite, s / J \ 1 elite, s / j \ 

(fc) and a. (fc): 



— elite, s 



(fc) 



3 

I -elite 

ap. 



/\ —elite / I T \ 

a)pj (k~ 1), 



(18) 



af *°'=(fc) = ad(fc)af *^(fc) + [1 - Qd(fc)] af'^{k - 1), (19) 

where a' is a smoothing constant parameter (0 < a' < 1) 
and Qd(fc) is a dynamic smoothing parameter at fcth itera- 
tion: 

Qd(fc) = a-Q(l-fc"')', (20) 

wit h < Q < 1 and q is an integer typica lly between 5 and 
10 (|Kroese. Porotskv fc Rubinsteinll2006l). 

Following ICaproni. Monteiro fc Abrahanil (|2009l l. we 
adopted a' — 1, a = 0.7 and q = 5 throughout this work. 
In addition, we also assumed A — 100, Aditc = 10 and 
fcmax ~ 800. It is important to emphasize that this value 
for fcmax was chosen so that, independently of the number of 
CE optimizations that are performed for a fixed Pprcc.obs and 
sense of precession, the remaining precession parameters will 
always converge to the same numerical value (variations not 
larger than a few per cent). For the present work, we run 
our code three times for each particular precession period 
and sense of precession. 



3.3 The merit function S 

The merit function ^(fc) transmits to the CE algorithm the 
best tentative set of precession parameters at the fcth iter- 
ation. Following Capro ni. Monteiro fc Abraham. (,2009 1. we 
have chosen ^(fc) as: 

S{k) = T(fc) + V { [Sc., (fc)]' + [Ss, (fc)]' + [Sr, (fc)]'} , (21) 



i = l 



where: 



(22) 
(23) 
(24) 



Sai{k) = Aa, - Aamodi(fc), 
^^^(fc) = AS^ - A5mcd,(fc), 
Sr,{k) = An - Ar-modi(fc), 

where Aa^ and A5i are, respectively, the right ascension and 
declination offsets of the jet knot i, Ar^ = Aa^ + AS^ . 

Our code generates, for each observation data fobs, pairs 
of right ascension and declination offsets based on different 
values of t^j , which are provided from a given jet-precession 
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model. The chosen model (RA,Dec)-pair is that which min- 
imizes the square distance between the precession helix and 
the observed (RA,Dec.) pair, also respecting the condition 
tobs.ej < ioba- Therefore, for each observed (RA,Dec.) pair 
there is always a tobs.oj that makes S^, + Sl, a minimum. 

As mentioned in ICaproni. Monteiro fc AbrahamI 

l|2009ll . the terms Sa^ and Sj. strongly constrain the 
instantaneous jet direction in the optimization process, 
while inclusion of Sr^ provides additional constraint on 
the modelling of the core-component distances as well as 
improving the convergence performance of the method. 

Jet kinematic modelling performed only on the plane of 
the sky has an extra potential difficulty: the jet components 
are not fully independent of each other, in the sense that 
unique identifications of jet knots also depend on estimates 
of their individual proper motions. However, identification 
problems related to fitting procedures, as well as observa- 
tions poorly sampled in time, may influence the follow-up 
of the components in time, which consequently might con- 
tribute to a misinterpretation of the data. Our CE modelling 
avoids such potential biases, purely analysing the time be- 
haviour of the sky position of the jet knots without any infor- 
mation concerning previous identification of components. To 
guarantee that our jet precession modelling will also predict 
apparent velocities similar to those observed in BL Lacertae, 
we included a penalty function T(k) in the CE o ptimiza- 
tion process l|Caproni. Monteiro &: AbrahamI [20091 ) . Based 
on independent estimates of the typica l apparent veloci- 
ties of the jet knots in BL Lacertae ( e.K., lMutel et al.|[l990l : 
iJorstad et al.ll200"5l : lLister et al.ll2009l '). we decided to include 
in the CE algorithm the following penalty function in equa- 
tion (21): 



T(fc) = 



10, 
0, 



if /3ir(fe)<2 
elsewhere, 



P^Tik) > 11, 



(25) 



where /3ob"(^) ''■iid I3'^^{k) are respectively the minimum 
and maximum values of the apparent jet speeds predicted 
by a given precession model at iteration k. It is important to 
emphasize that the choice of T(fc) = 10 mas^ is sufficient to 
guarantee that any tentative solution providing I3'^^{k) < 2 
or I3^^^{k) > 11 is statistically non-favoured during opti- 
mization. 

Note that the choice of 5* is based on the minimization 
of the quadratic distances between the observational data 
and those produced by the precession model, while T(A;) 
provides extra constraints on the parameters 7 and (j>o- 



4 THE JET PRECESSION MODEL FOR THE 
PARSEC-SCALE JET OF BL LACERTAE 

4.1 Observational data 

The pc-scale jet observational data of BL Lacertae anal- 
ysed in this work were gathered from the literature 
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Dec. offsets) were obtained at frequencies of 5.0, 8.4, 10.7 
15, 24 and 43 GHz, from 1980.93 to 2007.682, providing a 



Figure 1. The distribution of right ascension and declination 
offsets of the jet components of BL Lacertae. Grey circles rep- 
resent the whole data extracted from the literature, while the 
black circles are the data analysed by our precession model (core- 
component distances smaller than 2 mas). 



time monitoring of the jet activity in BL Lacertae of about 
26.7 yr. 

We show in Fig. [1] the corresponding spatial distribu- 
tion of the jet components of BL Lacertae on t he plane of the 
sky. As already noted in previous papers (e.g., Charlotiri99Gl : 
iDenn. Mutel fc Marscherll2000l : IStirling et al.ll2003l ). the pc- 
scale jet of BL Lacertae bends systematically to south-east 
after a core distance of about 3 or 4 mas, which definitely 
cannot be addressed exclusively by a jet precession phe- 
nomenon. Because of this, we decided to restrict our preces- 
sion analysis to distances below of 2 mas (see filled circles 
in Fig. [1}, avoiding the bending zone of the jet. 

Given the multiwavelength nature of the data used 
in this work, opacity effects are expected to affect the 
determination of the absolute core position, which in- 
troduce frequency-dependent shifts in the core-component 
dista n ces and proper motions (e.g., Blandford & Konig] 
I1979I : iLobanovl 1 19981 : iKovalev et al I I2OO8I I. In the case 
of a precessing jet, these corrections are time-dependent 
since the angle b etween the jet and the l ine-of-sight 
varies with time (ICaproni fc Abraham' '2 004al lbh. In the 
case of BL Lacertae, iMutel et al. (1990. ) found 0.3 mas 
for th e magnitude of the core- s hift b etween 5.0 and 10.6 
GHz. iDenn. Mutel fc Marschej (|2000l l estimated a core- 
compone nt shift between 22 and 43 GHz smaller than 
0.1 mas. ICroke fc Gabuzdal (|2008l ) calculated a core shift 
of 0.20 mas comparing 5.1 - and 7.9-GHz maps, while 



lO 'Sullivan fc Gabuzdal (|2009l ). using images at frequencies 
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Figure 2. Results for clockwise sense of precession (t = 1). Dependence of the CE optimized precession model parameters in terms of 
the observed precession period, as well as the value of the merit function (lower rightmost panel). Steps of 1 yr in precession period 
were considered, except for the interval between 11 and 14 years, during which differences between consecutive precession periods of 0.1 
yr were used. We also included in our analyses the value of 2.29 yr obtained by Stirling ct al. (2003). The merit function reaches its 
minimum at a precession period between 11 and 13.5 yr (dotted vertical lines). The solid and dotted vertical lines mark respectively the 
weighted mean and standard deviation of the precession period related to the merit-function's minimum (12.11 ± 0.65 yr). 



between 4.6 to 43.1 GHz, found core-shifts values below of 
0.43 mas. Considering these works, the mean value of the 
core-shift in BL Lacertae is roughly 0.2 mas, which will not 
influence significantly our jet precession modelling given the 
intrinsic uncertainties of the data, as well as the high frac- 
tion of 15-GHz data in our sampl^U- Because of that, we 
will hereafter neglect core opacity effect in our jet preces- 
sion analyses. 



4.2 Estimating the precession period and the 
sense of precession of the BL Lacertae jet 

The precession period and the sense of precession are quanti- 
ties not automatically optimized by our CE method. Thus, 
some extra procedu re is necessary to determine those pa - 
rameters. Following ICaproni. Monteiro k, AbraharnI l|2009l ). 



^ We run additional CE optimizations considering two subsets 
of the original data to verify quantitatively this claim. The two 
subsets correspond to 15-GHz and 15-1-43 GHz data (about 44 
and 83 per cent of the original data set, respectively). For 15-GHz 
data, typical changes in the values of the precession parameters 
in relation to our best model (see Table 1 in Section 4.4) are 
usually less than ~7 per cent (not exceeding ~25 per cent). For 
15-1-43 GHz data, the typical differences are less than ~3 per cent 
(~15 per cent in the worst case). The decrease in those values is 
a consequence of the increase of the time coverage of 15-1-43 GHz 
data (from ~ 8.3 yr to ~ 13 yr). 



we mapped the dependence of the jet-precession model fit- 
ting with the precession period and the sense of precession. 
We varied the value of the precession period measured in the 
observer's reference frame from 2 to 30 yr in steps of 1 yr 
for both clockwise and counterclockwise senses of precession, 
covering those value s previousl^y sugg ested in the literature 
(|Stirling et al.l [20031 : iTat evama 2009). An extra refinement 
step of 0.1 yr was adopted between 11 and 14 yr, in which 
the merit function presented a minimum, as will be discussed 
later. 

For each tentative precession period, our CE method 
searched for the best set of precession model parameters 
in the ranges: 0.95 ^ /3 ^ 0.999 (3.2 ^ 7 ^ 22.4), 
-180° ^ 770 ^ -150°, 0?1 ^ (/-o ^ 40°, 0° sC ifio sC 30° 
and ^ ro,s 5; 1. It is important to emphasize that 
these parameter ranges contain the jet paramet er values for 
BL Lacertae derived in previous works (e.g.. |Mutel et al.l 
1 1990"; 'Penn. Mutel fc Marsclieil l200d : [Stirling et all l2003l : 
Tatcvama 20o|)~ Note also that the hypothesis of a non- 
precessing jet is incorporated in our analysis, since ifio = 0° 
is one of the possibilities to fit the observational data. 

We show in Figs. [2] and [3] the optimized precession 
model parameters and the merit function versus the pre- 
cession period in the observer's reference frame, for clock- 
wise and counterclockwise senses of precession, respectively. 
Among all precession model parameters, 770 exhibits the 
smallest variation in terms of Pproc.obs (only about 4°). Ex- 
treme values for 7 and cj>o were found for some tentative 
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Figure 3. Results for counterclockwise sense of precession (t = —1). Dependence of the CE optimized precession model parameters in 
terms of the observed precession period, as well as the value of the merit function (lower rightmost panel). Steps of 1 yr in precession 
period were considered, except for the interval between 11 and 14 years, during which differences between consecutive precession periods 
of 0.1 yr were used. We also included in our analyses the value of 2.29 yr obtained by Stirling et al. (2003). The merit function reaches 
its minimum at a precession period between 11 and 13.5 yr (dotted vertical lines). The solid and dotted vertical lines mark respectively 
the weighted mean and standard deviation of the precession period related to the merit-function's minimum (12.11 ± 0.65 yr). 



precession periods; they correspond to 5 < Pprcc.obs < 14 
for clockwise precession, and -Pprcc,obs 5 or -Pprcc.obs ^ 19 
for counterclockwise precession. 

A well-defined minimum for the merit function is found 
at a precession period of about 12.1 years, independently 
of the sense of precession. The values of both minima are 
practically identical {S = 4.511 ± 0.109 mas^ and S = 
4.606 ± 0.072 mas^ for clockwise and counterclockwise pre- 
cession, respectively), making impossible to obtain informa- 
tion about the sense of precession from the merit function 
alone. The difficulty in det ecting the cor r ect se nse of pre- 
cession was also found by IStirhng et all (|2003l ). As men- 
tioned by them, the problem probably resides on the lack 
of observational data with larger core-component distances 
employed in the modelling. However, other constrains on 
the optimized precession parameters can be used instead, 
namely the interval of observed superluminal velocities, the 
boosting parameter, which depends on the Doppler factor, 
and the jet-to-counterjet flux density ratio H, which should 
be large enough to guarantee that the counterjet is not ob- 
served. This last quantity can be calculated from: 



1 -f /3 cos 0(rs) 
1 - 13 cos (P{Ts) 



p+a 



(26) 



where p is equal to 2 or 3 for a continuou s or clumpy jet, re- 
spectively (e.g.. lLind fc Blan dford''l985'). and a is the spec- 
tral index of the flux density distribution in terms of the 
observed frequency f {Sv oc «^ 



We plotted the dependency of the interval of predicted 
superluminal velocities, Doppler factor and jet-counterjet 
ratio as a function of the precession period in the ob- 
server's reference frame in Figs. U (clockwise sense of pre- 
cession) and [5] (counterclockwise sense). Looking carefully 
at PprGc,obs ~ 12.1 yr we can see that the clockwise jet- 
precession model predicts too low value of H (~ 300), which 
implies the possibility of counterjet detection in interfer- 
ometric radio images, given their typical dynamic ranges. 
However, such detection has never been achieved, which in- 
dicates that a clockwise jet-precession scenario is not ap- 
propriated for BL Lacertae. On the other hand, the coun- 
terclockwise jet-precession model predicts H ~ 3.5 x 10^, in 
agreement with available observational data. 

In addition, our counterclockwise precession model pre- 
dicts Doppler-boosting factor ranging from 8.8 and 9.4 (in 
contrast to the clockwise model that predicts a variation be- 
tween 1.2 and 1.7), in good agreement with independent es- 
timates of the Doppler-boosting parameter fo r BL Lacertae 
found in the literature that suggest ^ ^, 7 (e .g., |Jorstad et al.l 
I2OO5I : iHovatta et al.|[2009l : IWu et al]|201ll ). Therefore, the 
counterclockwise 12.1-yr jet precession scenario is favoured 
in the case of BL Lacertae. 

The 12.1-precession period found in this present work 
is not in agreement with the preces sion periods of 2.29 



and 26.0 y ears reported previously bv lStirling et al.l (|2003|) 
and Tatevama (2009), respectively. However, Mutel Sz DennI 
analysing a different set of 43-GHz images (includ- 
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Figure 4. Results for clockwise sense of precession (t = 1). Ap- 
parent velocity, Doppler-boosting factor and jet-counterjet ratio 
(assuming p = 2 and a = 0.8) predicted by our CE optimized 
model as a function of the precession period. Full circles represent 
those quantities calculated at (p = (po, while the dotted lines show 
their respective upper and lower limits allowed by the precession 
model. The solid vertical line marks the favoured precession pe- 
riod of 12.1 yr. 
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Figure 5. Results for counterclockwise sense of precession 
(t = —1). Apparent velocity, Doppler-boosting factor and jet- 
counterjet ratio (assuming p = 2 and a = 0.8) predicted by 
our CE optimized model as a function of the precession period. 
Full circles represent those quantities calculated at ip = ipg, while 
the dotted lines show their respective upper and lower limits al- 
lowed by the precession model. The solid vertical line marks the 
favoured precession period of 12.1 yr. 



ing some overlapping epochs with IStirling at allbooj ). con- 
cluded that changes in the structural position anglso of the 
radio core of BL Lacertae could present a periodicity of 12.1 
yr, as statistically significant as that of 2.29 yr. Note that all 
these previous works analysed a substantially shorter inter- 
val of observations (~3.1 years in .Stirling ct al.. .2003i , ~5.1 

^ Angle defined by the relative orientation between components 
CI and C2 detected in the 43-GHz VLBI images of BL Lacertae 
(see IStirling et "al]|2003h . 



years in lMutel fc Dennll2005l and - 12 yr in lTatevamall2009l ) 
in comparison to this present work (~26.7 yr), which might 
be responsible for such discrepancies. Another possibility is 
that the 2.29-yr period is related to some extra phenomenon, 
such as nodding modulation of the precession angle, for in- 
stance (see Section 6 for a more detailed discussion about 
this possibility). 
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4.3 The non-precessing jet scenario 

As mentioned before, tiie non-precessing jet liypotliesis 
[ifio = 0°) was also diecked by our CE teclinique dur- 
ing eacli op timization processes. Even tliougli tiie vali- 
dation tests l|Caproni. Monteiro fc Abraliamll2009l ) showed 
that the technique is able to determine whether preces- 
sion is taking place, we decided to perform an indepen- 
dent model optimization, forcing (po = 0°, and rerun the 
code to optimize the remaining jet model parameters. We 
found that the value of merit function (defined by equa- 
tion 21) in this situation corresponds to ~7.8 mas^, at 
least 1.5 times larger than that found for Pprocobs ~ 12.1 
yr. To estimate the statistical significance of such dif- 
ference, we applied the non-parametric two-sample two- 
dimensional K olm ogorov— Smirnov test, (2DKS test here- 
after: see e.g.. [ l^ea cock 1983; Fasano & Franccschini 198^; 
iPress et al.lll997l ). 

In brief, this test compares two-dimensional data sets in 
terms of the maximum cumulative difference between them, 
measuring the probability of these two data sets being drawn 
from the same parent distribution. In this work, the data sets 
are the observed and model predicted jet-knot position off- 
sets from the core, so that comparison is made using pairs of 
(Aoobsi, A(5obsi) ™<i (Aaniodj, A^modj), the ith data point 
of the respective data sets. Note that these two samples 
are not strictly independent since the model coordinates are 
calculated from a model that was constrained, on the other 
hand, from the observed positions of the jet knots. Although 
it is in contradiction to the underlying hypothesis of inde - 
pendence of the two analysed data sets, IWild et al.l (|2005l ) 
argue if the number of the observed data points is much 
greater than the number of the model parameters, such ef- 
fect should be small enough to allow the application of the 
2DKS test. This requirement is fulfilled in our data sets. 

Unfortunately, another potential drawback concerning 
the multidimensional generalization of the Kolmogorov- 
Smirnov test is related to the calculation of the max- 
imum cumulative difference in two or more dimen- 
sions, which is not uniquely defined (e.g . , Peacockl 1 19831 : 
iFasano fc Franceschinl 119871 : iPress et all 119971 ) It intro- 
duces a possible dependency in the determination of the 
maximum cumulative difference with the shape of the 
parent two-dimensiona l probability distribution. However, 
iFasano fc Franceschinil (| 19871 ) showed from Monte Carlo 
simulations that their parametrisation for 2DKS test is 
distribution-free in good approximation. In any case, the 
resulting probabilities obtained from 2DKS tests are not as 
rigorously reliable as in the one-dim ensional case, so that 
they must be interpreted warily (e.g.. |Press et al.]|l997l ). 

The application of the 2DKS test for our counterclock- 
wise 12.1-yr jet precession model leads to a proba bility 
of ~15 per cent. As suggested in IPress et al.l (|l997l ) and 
IWild et all (120051 ). a 2DKS probability above ~ 10-20 per 
cent can be used to claim reasonable agreement between the 
two analysed data sets. Conversely, this probability drops to 
0.028 per cent for the non-precessing jet scenario. Thus, the 
2DKS tests indicate that the probability of the observed 
and model data sets being drawn from the same parent dis- 
tribution is 15 per cent for our best precession model, but 
less than 0.1 per cent for the non-precession scenario. This 
demonstrates that our model provides a significantly bet- 



Table 1. The best precession model parameters optimized by our 
CE algorithm for a counterclockwise sense of precession (t = — 1). 
The uncertainties in each parameter correspond to la level. 



-Pprec,obs (y^) 


12.11 ± 0.65 


-Pprcc,s (yi") 


550 ± 247 




0.9824 ± 0.0087 


7 


5.35 ± 1.31 


VO (deg) 


-165.86 ± 0.16 


00 (deg) 


4.43 ± 2.16 


'PO (deg) 


0.51 ± 0.24 


"ro.s 


0.242 ± 0.008 



ter description of the collected observations for BL Lacertae 
than a non-precessing model, although some discrepancies 
between our best model and observational data still remain. 

4.4 Model parameters for 12.1-years 

counterclockwise precession period 

As presented in last section, the precession of the pc-scale 
jet of BL Lacertae occurs in a counterclockwise sense, re- 
specting a periodicity of about 12.1 yr in the observer's ref- 
erence frame. The CE optimized precession model parame- 
ters related to this period are given in Table [1] Precession- 
model parameters and their respective uncertainties were es- 
timated foUowing lCaproni et al.l (|201ll ) (equations 10 and 11 
in their paper). Note that the precession rate in the source's 
reference frame is substantially slower, as expected from rel- 
ativistic effects (see equation 11). 

Comparison between snapshots of the precession helix 
generated from the parameters listed in Table [T] and the 
right ascension and declination offsets of the observed jet 
knots of BL Lacertae can be visualized in the section enti- 
tled Supporting Information, published in the online version 
of this work. The distribution of the residuals in right as- 
cension and declination directions (generated from the dif- 
ference between the model-predicted and observed positions 
of jet knots) is compatible with zero-mean value in both di- 
rections. Their associated root mean square values (rms) are 
~ 0.12 and ~ 0.03 mas in right ascension and declination 
directions, respectively. It strongly indicates that the mean 
jet position angle predicted from our best model is quite well 
determined. 

Applying the same analysis for a non-precessing jet, it 
is possible to note an increase of the spread of the residuals 
in relation to that found from our best precession model. 
It is quantitatively corroborated by the slight increase of 
the rms in both directions (~ 0.15 mas in right ascension 
and ~ 0.04 mas in declination), reinforcing our conclusion 
that the non-precessing jet scenario is not suitable for BL 
Lacertae. I n relation to the pre viousl y published prece ssion 
models bv [Stirling et all (|2003l ) and iTatevainal (|2009l ). the 
amplitude of residuals is even higher in comparison with a 
non-precessing jet scenario and our best precession model. 
Indeed, this was alrea dy expected given t heir larger values 
of S (^ 17 ma s^ for IStirling et all |2003| and 15 mas^ for 
iTatevamal |2009| . much higher than our own value for no- 
precession case)If| 



In order to compare the merit function of our best jet- precession 
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According to our precession model, the expected 
apparent velocities of the jet of BL Lacertae are al- 
ways superluminal, ranging from about 3.4 to 4.0c. 
Components having suc h velocities have b een observed 
in previous works (e.g., Mutel et al.l 1990l: Jorstad et al.l 



I2OO5I : iGabuzda fc Cawthornell2003l : iLister et alj|2009ll , even 
though some quasi-stationary components (/3oba << 1), as 
well as f aster knots have been also detected in BL Lacer- 
tae (e.g.. I Jorstad et al.|[2005l : iLister et"al] 120091 '). It should 
not be interpreted as a serious problem of our model since 
some jet components of BL Lacertae exhibit non-ballistic 
trajectories, which is not taken into account in our present 
parametrization. Indeed, the theory of linear perturbations, 
as well as numerical simulations, has shown that jet preces- 
sion can_hiduce instabilities in the jet flow (Hardee 2000:, 
I200ll . l2002l ). The observational counterparts of these jet in- 
stabilities may be seen as non-ballistic knots in VLBI im- 
ages. Therefore, the scenario adopted in this work (ballistic 
motions plus jet precession) does not exclude the existence 
of helical motions in the jet of BL Lacertae. Note also that 
our predi cted apparent veloc ities are different from those 
found by IStirling et all (|2003l ) (/3obs ~ 6.7) and iTatevamal 
(|2009l ) (7 < /3obs < 10). 

Because of precession, jet viewing angle varies approxi- 
mately between 4° and 5° , which is in reasonable agreement 
with the values derived from the characteristics of the vari- 
ability in the multiwavelength light curv e s of BL Lacertae 
jHovatta et al.l [20091: iRaiteri et al] I2OI0I : ISavolainen et all 



12010!; 



Wu et al. 



20111 ). 



The predicted jet viewing angles 
from our precession model are between those found by 
IStirhng et all (|2003l ') (6° < < 12°) and iTatevainal (|2009i ) 
(1.5° <<^<3°). 

The jet inlet position angle on the plane of the sky 
ranges from -172° to -159°, providing a total variation of 
^ 13° which is app roxi mately half of tho se obtained by 
IStirhng et all (|2003f ) and iTatevamal (|2009l ). This angular 
amplitude means an amplitude variation of 0.46 mas at 2 
mas from the core, which is about twelve and four times 
greater than typical observational uncertainties in the right 
ascension and declination positions of the jet components, 
respectively. 

As mentioned previously, the Doppler-boosting parame- 
ter varies roughly from 8.8 to 9.4 in our model, which means 
that the jet-counterjet flux ratio ranges from 3.3x10^ to 
3.9x10^ (consi dering a continuous jet with a spectral in- 
dex of 0.8; e.g.. lLind fc Blandfordlll985l ). implying m a non- 
detected counterjet at pc s cales. T his result also supports 
the hypothesis claimed bv . Stir ling et al., (i2003) that compo- 
nent CI is not the manifestation of the counterjet. 

The observational data used to constrain our optimized 
precession model covers 26.7 yr of monitoring, which means 
an interval of about 2.2 times longer than the estimated 



mod el quantita t ively with those derived bv I Stirling et al" 
and ITatevamal l l2009l 'l.we generated precession helices based on 
the values of the parameters determined by those authors and 
compared them with the same observational data used in our 
work. Thus, these comparisons respect the same conditions used 
in this work (e.g. penalty function, jet region, frequency range, 
etc.), as well as their related small 2DKS test probabilities (0.039 
per cent for [ Stirling et al.l [20031 and 1.6 X 10~^ per cent for 
lTatevamall2009l) . 



precession period of 12.1 yr. At this point it is important 
to emphasize that our CE model technique is able to deal 
properly with short time s amplings, as can be seen in the 
validation tests published in lCaproni. Monteiro fc AbrahamI 
(2009). One of these tests employed synthetic data formed 
by no more than 50 sky positions (about one sixth of the 
data analysed in this work) and covered only 1.5 precession 
periods. Therefore, this strongly suggests that our results 
are not influenced by the relatively short time covering of 
the data in terms of the 12.1-yr precession period. 



5 THE UNDERLYING JET OF BL LACERTAE 

The jet precession period of 12.1 yr has no counterpart de- 
tected in analyses of the historical light curves of BL Lac- 
ertae conducted up to now. Quasi-periodic variations in the 
radio light curves of BL Lacertae have been reported in the 
literature, ranging roughly from 0.7 to 8 yr depe nding on the 
frequency and time cover age of the observations (|Kellv et al.l 
I2OO3I : IViUat"a et al.ll20o3 ). In optical domains the situation 
is almost the same: there is some eviden ce of a periodicity of 
about 7-8 yr but this is sti ll inconclusiv e jHagen- Thorn et al.l 
1 19971 : IVillata et all |2004| ). Similarly, iFan et all j 19981 ) re- 
ported a strong signature of a periodicity of 13.97 yr in the 
B-band light curve, a s well as shorter p eriods of 0.66 and 
0.88 yr, also found bv lWebb etahl (|l988l '). 

Although some works have shown the viability of jet 
precession in producing det ectable signatures in the li ght 
curves of some sources (^^e.g.. ICaproni fc Abrahamll2004al g). 



it should not be expected f or all precessing objects. As men- 
tioned bv iMutel fc DennI (120051 ). the main contribution to 
the observed flux of BL Lacertae may be due to the ex- 
tended jet, which is not varying its orientation and conse- 
quently may be masking any precession signature. Another 
possibility that cannot be excluded is that the underlying jet 
is intrinsically weak, so that even Doppler-boosting effects 
modulated by precession are not sufBciently strong to make 
its contribution detectable. 

We compared the historical 8.4-GHz and B bancj3 light 
curves with upper limits for the observed flux density of 
the underlying jel|f| predicted from our best model. We as- 
sumed simplistically that the intrinsic flux density of the 
underlying jet Ss did not vary along the observation time. 



fo llowing the relation 6'obs(tob s) = Ss[5(to 



(similarly 



to lCaproni fc Abrahamll2004bl). We assumed a ~ 0.0 for the 
radio regime llDenn. Mutel fc Marscheil I2OO0I ) and a = 2.0 
for the optical band ( Brown et al. 19891 ). To keep the under- 
lying jet contribution always below the measured total flux 
densities, Ss must be lower than 18 mjy and 0.30 fiJy at 
radio and optical frequencies, respectively. 

The small variation in the value of the Doppler-boosting 
factor (between ~ 8.8 and 9.4) implies a small variation in 



* Based on data taken and assembled by the WEBT collaboration 
and stored in the WEBT archive at the Osservatorio Astronomico 
di Torino - INAF (http://www.oato.inaf.it/blazars/wcbt/). 
^ The term underlying jet refers to a continuous or quasi- 
continuous distribution of plasma elements in the jet in which 
shocks and/or plasm ons (jet components) pr opagates (in a simi- 
lar sense as given bv lLind fc Blandfordlll985^ . 
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the amplitude of the observed flux density of the underly- 
ing jet (between ~ 0.18 and 0.23 mjy in the B band and 
from about 1.39 and 1.59 Jy at 8.4 GHz). We speculate that 
it might be the reason for the non-detection of such 12.1 yr 
periodicity from the previous statistical analyses of the light 
curves of BL Lacertae. Notwithstanding, the underlying jet 
contribution might be responsible for a plateau-like struc- 
ture (with mean values of about 0.21 mJy in the B band 
and 1.49 Jy at 8.4 GHz) underneath of the main flux vari- 
ations of BL Lacertae, which are probably produced by the 
jet components themselves; this c an be interpreted as either 
shocks moving down the jet (e.g. iBlandford fc Koniglll 19791 : 
iHughes. Aller fc AlleJ 19851. 19891) or plasrnons e lected from 
the core (e.g.. |Paulinv-Toth fc Kellermannll 19661 ). 



6 



JET PRECESSION AND NODDING 
MOTIONS IN BL LACERTAE 



[Stirling et al.l (|2003l ) analysed the structural position angle 
defined by the relative orientation between components CI 
and C2, as well as the polarization angle in radio and optical 
observations, finding a periodic modulation of 2.29 yr in 
both quantities. They interpreted those results as due to the 
precession of the jet inlet region. We explore in this section 
an alternative interpretation for this short period in terms of 
a nutation motion produced in a supermassive binary black 
hole system. 

iKatz et all (Il982l ) established the relation between the 
characteristics of a precessing accretion disc in close binary 
systems and the frequency ujn of short-term nodding mo- 
tions: 



(27) 



where cjorb is the orbital angular velocity of the secondary 
black hole. Note that Us should be negative since the induced 
precession is retrograde, in the sens e of being contrar y to the 
rotation of the accretion disk (eg.. lKatz et al.lll982l ). 

Assuming that the nodding oscillation has a period of 
2.29 yr in the observer's reference frame, we can use equa- 
tion (27) to calculate the value of the orbital period of the 
secondary black hole Porb.obs to produce the inferred 12.1- 
yr precession period. It implies Porb.obs = 7.4 ± 1.8 yr or 
forb.s ~ 335 ± 173 yr in the source's reference frame, the re- 
spective uncertainties having been obtained from error prop- 
agation of the involved periods. 

From the Porb.s, we can derive the distance between the 
two black holes Rps from the Kepler's third law: 



Rr, 



GMtot 2 

47^2 -^orb.s, 



(28) 



where G is the gravitational constant, and Mtot = Mp + Ms, 
the sum of the masses of the primary and secondary black 
holes, respectively. 

Several works have attempted to estimate 
the mass of the supermassive black hole in BL 



Lacertae |Fan. Xie & Bacon 




119991: IWu fc Urryl 


2OO2I: Marconi fc HuntI 2003 




Ghisellini et al. 20ld: 


CaDetti, Raiteri & Buttigliond 201C 


), which seems to have 



about 1.5 — 6 X 10 Mq. Assuming that Mtot corresponds 
to the mean value of those estimates, i.e Mtot = 3.75 x 10* 
we calculated the distance between the black holes 



from equation (28), providing Rps = 0.17 ± 0.07 pc, which 
corresponds to ~ 0.13 mas at the distance of BL Lacertae. 

To put a lower limit for the ratio between the masses 
of the secondary and primary black holes, q = Mg/Mp, it 
is necessary to verify in which conditions the orbital period 
of the secondary can produce the inferred accretion disc/jet 
precession rate. Assuming that jet precession of BL Lacer- 
tae is induced by torques in the primary accret ion disc due 
to a non-coplanar seco ndary bla ck hole (e.g., Katj 19971: 
Papaloizou fc Terauem|[l9 95: Lar wood)ll997l: Romero et al.l 



2OOOI : ICaproni fc AbrahamI l2004allbl: ICaprorii et al.l 1200^ 



and taking Tiprec ^ Rout ( Romero et al.ll2000l ). where Rpr, 
and -Rout are respectively the outer radii of the precess- 
ing part of the dis c and the disc itself, we can write 
ICaproni et al.ll2006l ): 



-fprcc.s \ ....^4 

— cos (fio^ 7; 

^orb.s / J 



\7 -2nJ 



+ 



1/3 



-1 3/2 



0.88g2/3/(q) 



(29) 



where n is the polytropic index of the gas (e.g., n = 3/2 for 
a non-relativistic gas and n = 3 for the rclativistic case), 
and the function f{q) has the form ([Eggleton. , 1983): 



0.49g 



2/3 



0.6.32/3 -Hln(l-f-Qi/3 



(30) 



The value of the left-hand of equation (29) is completely 
defined by our precession and nutation models, while the 
right-hand of the same equation depends only on the pa- 
rameter q. Assuming n — 3/2, equation (29) is satisfied 
only if g > 5.75, which implies Mp < 5.5 x 10^ Mq and 
Ms > 3.2 X 10* Mq. It is important to emphasize that there 
is no problem concerning secondary black hole to be more 
massive than primary one. The word primary in our con- 
text only means that the observed jet is generated from the 
primary black hole. It is not new at all, since some previ- 
ous papers had already claimed systems with more massive 
secondary black hole tha n primary ones (e.g.. lRomero et"al] 
|2000| : iBritzen et ai]|200ll ) . 

The time stability of such system can be verified by 
calculating the timescale raw for losses due to gravitational 
radiation from the expression |Begelman. Blandford fc ReesI 
ll980l : IShapiro fc Teukolskvll 19831 ): 



5^ 



raw 



256 (GMtot )3 



(31) 



M 



0' 



which is minimized when q = 1. 

For the supermassive binary black hole system consid- 
ered in this work, we have tqw ^ 70 Gyr, which is much 
longer than the age of the Universe. This indicates no signif- 
icant changes in the orbit of the secondary and consequently 
to jet precession rate during the interval of the observations 
used in this work. 

Concerning the intrinsic amplitude of the nutati o n mo - 
tion, our calculation based on equation (4) in iKat j (|l997l ) 
leads to a value of about 0?10. Its inclusion into our jet- 
precession model generates an additional amplitude varia- 
tion at the position angle of the jet on the plane of sky of 
about ±2? 5, which may be detected in the millimetre data, 
in principle. 

We can see that a putative supermassive binary black 
hole system with the physical characteristics described 
above can drive not only the jet-precession period found in 
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this work, but also a nutation with a period of 2.29 yr, which 
might be responsible for variations in the structural position 
angle of the pc-scal e radio core, a s well a s in t he polarization 
angles reported by [Stirling et al.l (|2003l ) and lMutel fc DennI 
(|2005l ). 



frame is approximately 335 yr, which means a distance be- 
tween primary and secondary black holes of about 0.17 pc. 

Further monitoring of the pc-scale activity of BL Lac- 
ertae is necessary to confirm the validity of our precession 
model, as well as the supermassive binary black hole scenario 
proposed in this work. 



7 CONCLUSIONS 



[Stirling et al] (|2003h and lTatevamal (|2009l ') have claimed the 
existence of precession motions in the pc-scale jet of BL 
Lacertae of 2.29 and 26 yr, respectively. In this work we 
revisited this issue, investigating temporal changes of the 
observed right ascension and declination offsets of the jet 
knots in terms of our relativistic jet-precession model. 

Our precession model is characterized by seven pa- 
rameters: precession period, jet bulk velocity, aperture an- 
gle of the precession cone, the angle between the cone 
axis and the line of sight, position angle of the cone axis 
on the plane of the sky, sense of precession, and preces- 
sion angular phase. These parameters were optimized via 
heuristic cross-entropy (CE) method, comparing the pro- 
jected precession helix with the two-dimensional position 
of the jet components on the plane of the sky, follow- 
ing ICaproni. Monteiro fc AbrahamI (|2009l ) . The search for 
the best precession model for BL Lacertae is performed 
considering parameter ranges that contain th e jet param- 
eter values derived in previous wor k s fe.g.. | Mutel et aU 
1990l: iDenn. Mutel fc MarscheJ I2OO0I : [Stirling et all l2003l : 



Tatevarnaf 20091 ). as well as the case of a non-precessing jet 



{^0 = 0°). 

Our optimized best model is compatible with a jet hav- 
ing a bulk velocity of 0.9824c, which is precessing in a coun- 
terclockwise sense with a period of about 12.1 yr at the 
observer's reference frame 550 yr at the source's refer- 
ential framework), and changing its orientation in relation 
to the line of sight between approximately 3?9 and 4?9. The 
position angle of the precession cone axis on the plane of 
the sky is about -166°. It is important to emphasize that 
these precession parameters lead to 5 ~ 4.6 mas'^, which is 
significantly lower than the values obtained for both a non- 
precessing model and the previously published precession 
models. 

Searches for periodic variation in the historical light 
curves of BL Lacertae have not revealed any signature of 
flux variability occurring at a periodicity of 12.1 yr (e.g. , 
Webb et al]|l988l:lFan et al.ll 19981: iHagen-Thorn et al.lll997l : 



Kellv et al.ll2003l : IVillata et al.ll2004l '). We speculate that the 
non-detection of the periodicity of 12.1 yr might be at- 
tributed to small variations in the amplitude of the Doppler- 
boosting factor predicted by our model, which implies small 
variations of the observed flux density associated to the un- 
derlying jet. Nevertheless, the underlying jet contribution 
might be responsible for a plateau-like structure underneath 
of the main flux variations of BL Lacertae. 

Assuming that jet precession has its origin in a super- 
massive binary black hole system, we show that the 2.3- 
yr periodic variation in the structural position angle of the 
VLBI core of BL Lacertae reported by Stirling et al. is com- 
patible with a nutation phenomenon if g > 5.75. The orbital 
period of the secondary black hole at the source's reference 
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8 SUPPORTING INFORMATION 

Additional Supporting Information may be found in the on- 
line version of this article: 

This additional material shows the projections on the 
plane of the sky of the precession helices generated from the 
model parameters listed in Table [1] over the entire observa- 
tional period used in this work. We also display the 311 sky 
position of the jet knots of BL Lacertae that constrained our 
precession model. Note that precession helices describe the 
general trend of the observational data appropriately. 
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Figure 6. Changes in the jet orientation of BL Lacertae along the time due to the precession model parameters listed in Table ^ Solid 
lines represent snapshots of the precession helix, while the positions of the jet knots and their respective uncertainties are shown by black 
dots. Numbers at upper right corners of the panels are the corresponding observation epochs. 
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Figure 6. Continued. 
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Figure 6. Continued. 
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Figure 6. Continued. 
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